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Abstract 

We show how to compute chiral logarithms that take into account both the 0(0^) taste-symmetry 
breaking of staggered fermions and the fourth-root trick that produces one taste per flavor. The 
calculation starts from the Lee-Sharpe Lagrangian generalized to multiple flavors. An error in 
a previous treatment by one of us is explained and corrected. The one loop chiral logarithm 
corrections to the pion and kaon masses in the full (unquenched), partially quenched, and quenched 
cases are computed as examples. 

PACS numbers: 12.39.Fe, 11.30.Rd, 12.38. Gc 
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I. INTRODUCTION 



For simulating fully dynamical lattice QCD at light quark masses, staggered (Kogut- 
Susskind, KS) fermions have the advantage of being very fast relative to other available 
methods ^|. In addition, an exact chiral symmetry for massless quarks is retained at finite 
lattice spacing. However, the advantage in speed of KS fermions may be offset by systematic 
issues: on present realistic lattices {e.g., recent MILC simulations with a ^ 

0.13 fm), the KS taste^ violations are not negligible. Indeed, despite the fact that the MILC 
simulations use an improved ("Asqtad") action that reduces taste violations to (9(Q;|a^), 
these effects can still introduce significant lattice artifacts. 

Since one can control the taste of the external particles explicitly in the simulation, 
taste- violating artifacts show up primarily in loop diagrams. In particular, any quantity or 
computation that is sensitive to chiral (pseudoscalar meson) loops can be expected to show 
large artifacts at current lattice spacings. In order to perform controlled chiral extrapolations 
and extract physical results with small discretization errors from staggered simulations, it 
is necessary to include the effects of taste violations explicitly in the chiral perturbation 
theory (XPT) calculations to which the simulations are compared. The goal of this paper is 
to develop such a "staggered chiral perturbation theory" (SXPT). 

One can think of the MILC simulations as introducing flavor with separate KS fields for u, 
d and s quarks. The 4 tastes for each field are then reduced to 1 by taking the fourth root of 



the quark determinants for each flavor.^ The theory with vDet does not have a local lattice 
action, and there is some concern that non-universal behavior may thereby be introduced 
in the continuum limit. If we are able to show, by comparing simulations to SXPT forms, 
that the staggered theory produces the expected chiral behavior in the continuum limit with 



controlled 0{a?) errors, it should go a long way toward easing worries about the vDet trick. 
A starting point for any SXPT calculation is the work of Lee and Sharpe [G] , who derived 



^ We use the term "taste" to describe the staggered symmetry induced by doubhng; the taste symmetry 
becomes SU{4)l x S'C/(4)fl in the massless, continuum hmit, but is broken at 0{a?). We reserve the term 
"flavor" for true (u, d and s) flavor. 

^ Since m„ is always chosen equal to rud in the MILC simulations, one actually uses a slightly simpler 
procedure in practice. Only two KS fields are introduced, and the square root of the u, d determinant is 
taken. However, assuming algorithmic effects (step-size errors, autocorrelations) are under control, the 
two approaches are equivalent. We therefore prefer to consider the conceptually simpler case where each 
KS field represents a single flavor. 
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the 0{a'^) chiral Lagrangian for a single KS field (1 flavor, 4 tastes). In Ref. a general- 
ization of the Lee-Sharpe Lagrangian to multiple quark flavors was introduced to calculate 
chiral loop effects. However, there are subtleties in the generalization that were not appreci- 
ated in Ref. 0], leading to errors in the multi-flavor chiral Lagrangian and hence in the flnal 
chiral-logarithm formulas. These same subtleties also turn out to have implications even for 
the tree-level comparison (in Ref. jsl]) of the 1-flavor theory with simulations. 

Below we will follow the outlines of the three-step procedure introduced in Ref. [3], which 
we restate here for completeness: 

1. Generalize the Lee-Sharpe Lagrangian to correspond to n staggered quark flelds, re- 
sulting in a (broken) SU{4n)L x SU{An)ji chiral theory. Where convenient, we will 
specialize to the case of interest, n = 3. We call the n = 3 theory the "4 + 4 + 4" 
theory, since it has three flavors, each with four tastes; its symmetry is a broken 
SU{12)l X SU{12)ji. 

2. Calculate one loop quantities (such as m^^) in the 4+4+4 theory. 

3. Adjust, by hand, the result to a single taste per flavor in order to correspond to the 



physical case (and to simulation data). This adjustment corresponds to the yDet trick. 
It requires an understanding of the correspondence between the meson diagrams at 
the chiral level and the underlying quark diagrams and is basically the "quark flow" 
technique of Ref. Q]. For non-degenerate quark masses, we call the adjusted case 
the 1 + 1 + 1 theory; when we take m„ = rrid = mi (which corresponds to the MILC 
simulations) we call it the 2 + 1 theory. 

The difficulties in Ref. [3] arose in step 1. Fierz transformations were used to simplify 
the flavor structure in the taste-symmetry breaking potential. However, Ref. 0] had already 
employed Fierz transformations to simplify the form of this potential. The two transforma- 
tions turn out not to be compatible. In the Lee-Sharpe case, there was only one flavor, so 
this was not an issue. By properly taking into account the mixing of the flavor indices, we 
flnd that two of the six terms in the symmetry-breaking potential of Ref. |3] are incorrect. 

Another difference with Ref. is that there n was taken to be 2, and step 3 was modifled 
to adjust the ti, o? loops according to a v^Det, rather than a v^Det trick. This was due to the 
fact that Ref. |7[ took = nid from the beginning. However, the entire procedure is much 



clearer if every quark flavor is treated equivalently. Further, we will see that it is important 
to be able to treat directly charged pions [e.g., ud) that are composed of two independent 
flavors transforming under an exact lattice flavor symmetry (when m„ = m^^). Finally, the 
calculation is actually simpler when we keep all three quark masses unequal. The fact that 
the Goldstone charged pion mass squared must then have an overall factor of m„ + gives 
a very useful check on our calculation. 

Generalizing the taste-breaking potential properly has lead us to realize that flavor- 
neutral mesons in certain taste-nonsinglet channels can mix at tree-level due to "hairpin" 
diagrams. We can now see that such diagrams are present even in one-flavor XPT j^; 
their effects have however not been appreciated previously. The coefficients of the hairpin 
diagrams that arise here are new parameters in the chiral theory and have to be fit with 
simulation data or determined perturbatively. 

This paper mirrors the format of Ref. [3]. In Sec. m we generalize the Lagrangian of Lee 
and Sharpe, properly taking into account the flavor and taste structures involved. Sec. |in] 
discusses the calculation of the one loop chiral logarithms for the flavor-nonsinglet Goldstone 
meson mass in the 4+4+4 theory. It is convenient at this point to generalize the calculation 
to the partially quenched case, where the valence and sea quark masses are completely non- 
degenerate. The results are actually most simply expressed in this case, since there is a 
clear distinction between valence and sea quark effects, and no degeneracies arise that lead 
to cancellations. We then make the transition to the 1 + 1 + 1 theory in Sec. IIVI We write 
down results for both the partially quenched and "full" (equal valence and sea quark masses) 
cases, focusing primarily there on features which are different from Ref. ~\. The results for 
the quenched chiral logarithms are discussed in Sec. El Section |VT1 adds in the analytic terms 
and gives a compendium of flnal results, in full, partially quenched, and quenched cases. In 
the full iriu = TTLa = mi (2 + 1) case, the results from Sec. IVII have already been reported in 
Ref. We conclude with remarks about other uses for SXPT in Sec. IVIII An Appendix 
gives some additional details about the symmetries of the theory and briefly discusses the 
possible existence of a heretofore unknown phase of the staggered theory. This possibility is 
however apparently unrealized for physical values of the quark masses. 
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II. GENERALIZATION OF LEE-SHARPE LAGRANGIAN 



Lee and Sharpe p describe pseudo-Goldstone bosons with a non-linearly realized 
SU{A)l X SU{A)f{ symmetry, which originate from a single KS field. This KS field describes 
four continuum tastes of quarks. 

The 4x4 matrix S is defined by 

16 

S = exp(i0//), (f) = ^7raTa (1) 

a=l 

where the Ha are real, / is the tree-level pion decay constant (normalized here so that 
fjr ~ 131 MeV), and the Hermitian generators Ta are 

Ta = {^5, i^t,5, i^t^u, 6}- (2) 

Here we use the Euclidean gamma matrices with C^f^,^ = ^^^^ (/i < z/ in eq. (0)), ,^^5 = ^^^5, 
and ^/ = / is the 4x4 identity matrix. The field S transforms under SU{4)l x SU{A)ii as 

As discussed in Ref. 0, we will keep the singlet meson tt/ oc tr0 in this formalism. Due 
to the anomaly, the singlet receives a large contribution (which we will call itlq) to its mass, 
and thus does not play a dynamical role. Lee and Sharpe do not include this field in their 
formalism, which is equivalent to keeping the singlet in and taking mo — * 00 at the end of 
the calculation We keep the singlet here since in the generalized case of n KS fields, 
it is only the SU{4n) singlet that is heavy. In the mo — >■ 00 limit, the other SU{A) singlets 
will still play a dynamical role. 

The (Euclidean) Lee-Sharpe Lagrangian is then^ 

1 9m2 
^(4) = ytr(9^S9^St) - -/xm/2tr(S + St) + ^{^jf + a'V, (3) 

where /i is a constant with units of mass, and V is the KS-taste breaking potential. Correct 



^ Aside from the TOq term, we need not worry about wj dependence in this Lagrangian, since we are taking 
the niQ — > 00 hmit. It is only in the quenched case (Sec. 0), where we are unable to take the mg — > 00 
limit, that we will have to examine other tt/ terms. 
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through 0{a'^,m) in the dual expansion in and m, we have 



1 



k=l 



+C2^[tr(S2)-tr(e5Ee5S) + /i.c.] 
+C3^$^[tr(e.Se.S) + /i.c.] 

V 
V 
V 

+C6 ^tr(e^,Se,^St). (4) 

The 16 pions fall into 5 5*0(4) representations with tastes given by the generators Tq. This 
comes from the "accidental" 50 (4) symmetry of the potential V. We can determine the 
tree-level masses of the pions by expanding eq. Q to quadratic order: 

<=2/.m + ^5^,, + a^A«(CB), (5) 

where B e |5, /i5, /iz/(yU < z/), yU, /}. The A*^^^(^b) term comes from the V term, and is given^ 
in Refs. p, |7[ as: 



A^'^ (^m5) = (Ci + O2 + 3O3 + O4 - O5 + SOe) 
16 



A(^)(e5) = 
16 

16 

1 6 

A^'H^m) = (Ci + C2 + O3 + 3O4 + O5 + SOe) 
1 6 

A^'HO) = -^(403 + 404). (6) 
The vanishing of A(^''(^5) is due to the taste nonsinglet f/4(l) symmetry 

S ^ e^^«^Se^^«% (7) 
which is unbroken by the lattice regulator, making tts a true Goldstone boson. 



^ In Refs. 1^, Ldl, these corrections are denoted as A(^b). When we generahze to muhiple KS flavors, we 
wiU wish to distinguish this single flavor A'^^)(^b) from the n-flavor A(^b) 
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We now wish to generalize to the case of multiple KS fields. In Ref. for two KS quark 
fields, this was accomplished by promoting S and the mass matrix to 8 x 8 matrices. In 
the general case of n KS fields, which we discuss here, these become 4n x 4n matrices. The 
kinetic energy and mass terms are correctly given in Ref. The only difficulty arises in 
generalizing the taste-symmetry breaking potential (or equivalently the taste matrices ^b)- 
The generalization of V in Ref. 3] uses a Fierz transformation on the various four-quark 
operators to bring them into a "fiavor unmixed" form as follows: 



(8) 



where q is the quark field, z,j are SU{n) fiavor indices, 75 and 75/ are spin matrices, and 
and ^T' are taste matrices^. Treating the taste matrices as spurion fields, we see that for 
fiavor unmixed 4-quark operators, the ^ are singlets under the fiavor SU (n) symmetry. We 
can thus make the replacement: 



B 



iB 





y : ; ; 



(9) 



where the e}b^ are An x An matrices, and the E,b on the right hand side are still 4x4 taste 
matrices. 

Lee and Sharpe, however, already use Fierz transformations on the operators in Appendix 
A of Ref. jQ] to ensure that the final six operators in eq. (jD) are all single-trace objects. We 
now find that the transformation used in Ref. does not keep the operators in the same 
single-trace form. 

To see this, let us first assume we have made the replacement (jH)) in the taste-symmetry 
breaking potential. The operators O2 and are then not invariant under axial rotations of 
the individual fields. For example, consider a taste Ua{X) transformation on a single fiavor 



^ In Ref. |g, these are referred to as KS-flavor matrices and denoted by £,f and S,f' 
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only: 

/ 



^5 ■ ■ ■ 
1 ••• 
1 ■■■ 



(10) 



y : : : .J 

where S5 is a 4n x An matrix, shown here as composed of 4 x 4 blocks. It is simple to 
verify that the operators Oi, O3, C4, and are invariant under eq. (fTUj) . However, using 
gifi^s _ QQg^ _|_ sin 6', one finds that (^2 and are not invariant, and thus are not the 
correct generalization of the Lee-Sharpe terms to n flavors. 

One approach to generalizing the Lee-Sharpe Lagrangian correctly is therefore to consider 
all the different ways that the flavor indices on the various E fields in eq. (jl)) can contract. 
To do this, we write everything as 4 x 4 matrices and show the flavor indices explicitly. For 
example, the form of O2 from Ref. |^ can be written as: 

0^'^°^^'^'^' = ^[tr(E,,E,,) - trfeE^.^sS,,) + h.c], (11) 

where .^5 the 4x4 object, and i and j are the SU{n) flavor indices, to be summed over. 
Another SU{n) invariant we can create with this operator is: 

O2 = ^[tr(Ei,E,,) - tT{^,J:u^5^n) + h-c]. (12) 

One can easily see that this operator is invariant under eq. ()1U|). 

By starting with the other operators in eq. (^, we can similarly find other correctly 
generalized terms. This would for instance alter O5 along the same lines as eq. (fT^ . However, 
a problem with this approach is that it is difficult to ensure that the most general taste- 
violating potential is generated. For example, the operator tr(EjjE]^ — ^sEjj^sE]^) is invariant 
under eq. (fTUj) but is not easy to find starting with eq. (jH). That is because Lee and Sharpe 
have already set tr(EE''") = const, in arriving at their Oi. 

A more direct way to find the final form of the taste-breaking potential involves starting 
from the quark level and using the original analysis of Lee and Sharpe instead of their 
final result. At the quark level, gluon exchange can change taste and color, but not fiavor. 
Therefore the taste-violating 4-quark operators are composed of products of two bilinears. 
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each of which is a flavor-singlet, as in eq. (jH))- The 4-quark operators may be mixed or 
unmixed in color. ^ 

To 0(0^) in the dual a^,m expansion, the taste-breaking operators can be computed in 
the chiral limit. Since gluon emission does not change chirality, each bilinear is separately 
chirally invariant. The only such bilinears are vector and axial vector in the naive theory, 
which correspond to "odd" operators in the s tagg ered theory (operators in which quark and 
antiquark flelds are separated by 1 or 3 links) Thus only the odd-odd 4-quark operators 
in Appendix A of Ref. are relevant to us here. Each such operator can occur in color mixed 
and color unmixed form, but that does not affect the correspondence to SXPT operators.^ 
The even-even operators of Ref. were obtained by Fierzing the odd-odd operators and 
may be ignored: They correspond to flavor-mixed 4-quark operators. 

The above reasoning implies that the arguments in Ref. [7| were in fact correct, but only 
if the replacement eq. Q is implemented before the Fierz transformations in Ref. jsl that 
put the chiral operators in single-trace form. Writing the potential as V = W + W, we then 
obtain: 

= = CiTr(ef Ee^")st) 

k 

+C3^$^[Tr(e(")Ee(")S)+/..c.] 

V 

+C4^E[TKei?Se(?S)+/^.c.] 

V 

+^6 $^Tr(ei:)se;st) (13) 
-U'^Y.^,.0,. = C2yi$^[Tr(e(")S)Tr(e(")S) + /..c.] 

k' V 

+^24E[Tr(dfs)Tr(e(?S) + /^.c.] 

V 

+C5y^EtTr(e^^S)Tr(e(")st)] 

V 

+C54E[TK^i?S)Tr(e£^St)], (14) 



2 

V 



6 In Ref. color-mixed operators are Fierzed to put them in a standard, color-unmixed form. But this is 
precisely what we do not want to do here because it would mix the flavor indices. 

^ The color structure does affect the coefficients of the SXPT operators, but since these coefficients are 
arbitrary at the chiral level anyway, color mixing is irrelevant here. 
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where Tr is the full 4n x 4n trace, and t 



le (^^^ are An x An matrices as in eq. 



The terms 



that comprise U were found in Ref. j?]. Now, however, there are no terms that directly 
correspond to the operators O2 and O^. Instead, we have the four terms in U'.^ It turns out 
that only two combinations of the four constants in U ' enter in the 1-loop result: C2V ~ C'sy 
and C2A — C^A- The terms corresponding to C2V + Csy and C2A + C^a do not appear at 
this level. 

Note that the "accidental" SO {A) symmetry of the one-flavor theory jcj survives in 
eqs. (\Vi\ji and ()14|) . as seen by the fact that the the taste indices are contracted in a "Lorentz 
invariant" way. This implies that the degeneracies of the one-fiavor theory will also appear 
in the n-flavor case: all four taste-vector pions of a given flavor will be degenerate, as will 
all taste-tensors, etc. See the Appendix for further discussion. 

For n KS flavors, S = exp(i$//) is a 4n x An matrix, and $ is given by: 



^ U 71+ K+ 

71- D 

K- S 

V • • • 



(15) 



where U = Xla=i UaTa, etc., with the Ta from eq. Q. The component fields of the diagonal 
(flavor-neutral) elements {Ua, Da, etc.) are real; while the other (charged) fields are complex 
(tt^, K^, etc.), such that $ is Hermitian. Here the n = 3 portion of $ is shown explicitly. 
The mass matrix is now generalized to the An x An matrix 



M 



mj 
mJ 
msl 



\ 



\ ■ ■ ■ 

where again, the portion shown is for the n = 3 case. 
Thus, our (Euclidean) Lagrangian becomes: 

C = ^Tr(9^S9^St)-i/i/2Tr(A<S + ^St) 
o 4 



(16) 



2m, 



^Uj + Di + Si + ---)^ + a^V, (17) 



® The combination O2V + O2A can be Fierzed into the correct version of O2, eq. H12|l . and similarly for 
O5A — O5V Siiid the correct version of O5. The other linear combinations are new here, but could have 



been Fierzed into other operators of Ref. |3| if there were no flavor indices 
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where the rn^ term includes the n flavor-neutral fields and V = W + W is given in eqs. ()13|) 
and (HI). The ^^"^ in V are block-diagonal 4n x 4n matrices, as in eq. 0- 

When the masses vanish, the chiral Lagrangian, eq. (jl7p . has a flavor SU{n) vector 
symmetry and the individual Ua{^) symmetries for each flavor, both of which were used 
above, as well as overall fermion number conservation. These symmetries actually extend 
to a U{n)i X U{n)r "residual chiral group," although this full symmetry is not particularly 
important to us in the present context. Details are relegated to the Appendix. 

Expanding C to quadratic order in meson flelds, the potential U gives different masses 
to different taste mesons, but because it consists entirely of single-trace terms, the contri- 
bution is independent of the meson flavor. However, since U ' consists of two-trace terms, it 
contributes only to the masses of flavor-neutral mesons, and in particular only those with 
vector and axial vector tastes. Thus, even at tree- level and with rrtu = frid, a vr"*" of a 
given taste receives different mass corrections than a neutral U or of the same taste. In 
simulations, disconnected propagators for taste-nonsinglet pions (including the Goldstone 
pion) have invariably been dropped. This implies that simulations describe 7r+ mesons, not 
those constructed from a single flavor, which would have disconnected contributions. The 
comparison in Ref. |S| of the 1-flavor SXPT tree-level results to simulations is therefore not 
justifled, although almost all of the conclusions of Ref. survive a revised treatment. 

We thus want a chiral theory with both u and d quarks, even if we are interested in the 
rriu = md case. This is the primary reason that we consider the 4+4+4 theory here rather 
than the 4+4 theory of Ref. [3].^ 

From eq. (fTTjl . the tree-level masses of the mesons are: 

"^Ms = /^("^« + '^b) + a^A(^B), (18) 

^ We remark however that it would in principle be possible to extract the 7r+ results from a calculation 
in a (partially quenched) 4+4 theory. 



11 



where a and b refer to the two quarks which make up the meson M, and we have defined: 



A(e5) = Ap = A«(e5) = o 

16 



AiC^^) = AA = ^iCi + 3C3 + C4 + 3C6 



16 

A(e^,) = At = A(i)(e;..) = (2^3 + 2C4 + 4C6) 

AiQ = Av = ^iC^ + Cs + 3^4 + 3^5) 

1 6 

A(6) = A, = A(^)(6) = (4^3 + 4C4) . (19) 

Note that the mg terms and the terms from U ' are not included in these masses. Those 
terms, which affect only flavor-neutral mesons and give non-diagonal contributions in the 
basis of eq. (fT3j) . will be treated as vertices and summed to all orders below. Thus the 
corrections in eq. |T8|l are flavor independent. 

Simulations with the "Asqtad" action ^ give approximately equal splittings of the mass- 
squares of various taste mesons in the order M5, M^s, M^,^, M^, M/. From eq. (jTHl), this 
indicates that C4 is the dominant coefficient, a conclusion first noted in Ref. y. 

Upon expanding U ' in eq. (fT^ to quadratic order, we find a two-point vertex mixing the 
taste- vector, fiavor-neutral mesons (f/^, -D^, etc.): 

1 6 

-c?j^{C^v-C^v)^-aH'y . (20) 

In other words, there is a term +^^(f/^ + -D^ + "S"^! + ■ • 0^ in L. The vertex in the chiral 
theory is shown in Fig. QJa); while the corresponding underlying quark diagram is shown in 
Fig. ^b). There is also a vertex mixing the taste-axial, fiavor-neutrals (?7^5, -0^5, etc.): 

1 6 

-a2-^(C2A-C^5A) = -aX, (21) 

i.e., a term +^^(t/^5 + -D^^s + 5'^5 + ■ ■ in L. Similarly, the mg term in L produces a 
vertex — 4mo/3 between the taste-singlet, fiavor-neutrals (f//, Di, etc.). 

We thus have to resum the fiavor-neutral propagators in three cases: taste-vector, taste- 
axial, and taste-singlet. The methods of Appendix A in Ref. [3| allow us to calculate 
the full fiavor-neutral meson propagators easily and write them explicitly in terms of the 
true propagator poles (mass eigenstates). Here we sketch a few steps in this process. For 
concreteness we focus explicitly on the taste-vector case, although the taste-axial case is 
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obtained simply by replacing V with A in the equations below. The taste-singlet (^q) '^^^^ 

can be calculated similarly, although a more standard approach is also possible. We write 
the full inverse propagator as: 

Gy^ = Coy + (22) 

with 

{Goy)MN = (g^ + ml[^)6MN, (23) 

Hl,^ = a%, \fM,N. (24) 

Here and below we use V for generic taste- vector states, rather than the index /x. The indices 
M and N refer to the flavor- neutral mesons in the original basis of eq. ()15|). with tumv and 
mjsfy the "unmixed" masses from eq. (fTHj) {i.e., without including the mixing of eq. (j2ni))- 
For example, in the n = 3 case, these mesons are Uy, Dy, and Sy- Using Ref. we then 
find that: 

Gy = Goy + T)^ 

is the part of the taste- vector fiavor-neutral propagator that is disconnected at the quark 
level {i.e., Fig.[T]plus iterations of intermediate sea quark loops). We can write this explicitly 
in terms of the masses as: 



1 det(Gnv) 



rrir 



= ^=F • (27) 

(?' + m^M^W + ^''Ny)\{{<f + 

F 

Here L, like M and N , labels the unmixed fiavor-neutral mesons in the original basis 
are the poles of G^y); while F indexes the eigenvalues of the full mass matrix are poles 
of Gy). For n = 3, we name the corresponding full eigenstates in the taste- vector case Tiy,rjy, 
and rjy in analogy with the physical, fiavor-neutral, taste-singlet eigenstates. We emphasize, 
however, that all these taste-nonsinglet particles (including r]'y and the corresponding taste- 
axial particle rj'j^ are physically merely varieties of "pions:" pseudoscalars that do not couple 
to pure-glue states in the continuum limit, unlike the real r}' . 
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It is easy to generalize eqs. ()25l27p to incorporate partial quenching. Iterating Fig.mb) to 
determine the full propagator generates internal quark loops. Only sea (unquenched) quarks 
are therefore allowed in this iteration. Thus, if the number of sea quarks is n^, the product 
over L in the numerator of eq. (j27|) includes only the ns unmixed flavor- neutral mesons built 
from these sea quarks. Likewise, only the ns full eigenvalues are included in the denominator 
product over F. The external mesons M and A^, however, may be any flavor-neutral states, 
made from either sea quarks or valence quarks. Similarly, in the quenched case I^](^7v 
simply 

,T-«y,quench 2 r/ ^ 



(g2 + m|,^^)(g2 + m^^^ 



Below we will also need the relation 

det(G^i) _ _ 



tr(Goyi^^) 

(29) 



det(Go-y) 

Here the sum over L is again over the unmixed flavor-neutral mesons in the original basis. 
(In the partially quenched case, only mesons made from sea quarks are included in the 
sum.) This relation allows one to transform between the result ()27|1 and the form |13] one 
gets directly by iterating the 2-point vertex, eq. (jSOI)- 

Equations through apply explicitly to the taste-vector case; to get the taste- 
axial case, just let V ^ A. These formulas can also be used for the taste-singlet (/) channel 
with the replacement a^Sy —>■ Aml/S. We get: 

'i (p2 + ml,^ ) (p2 + m%^)[ {p' + m'p^ ) 

where L and F have the same meaning as in eq. (jlTj). The ttlq oo limit in the 4 + 4 + 4 
case is easily obtained, if desired, using m^^ = 4mQ for large ml. The rjj then decouples. 
However, we prefer not to take the ttLq oo limit at this stage, because the form of the 
result is then slightly different in the 4+4+4 and 1+1+1 cases, as we will discuss in Sec. IIVI 
In the quenched case, the product over sea quark states in the numerator and denominator 
of eq. (jHIH) are omitted. Of course, ml now cannot be taken to infinity, and the rj'j does not 
decouple. It is therefore necessary to consider possible additional rjj dependent terms in our 
Lagrangian. As discussed in Refs. j?! and 14 1, one can do this simply by making the the 



14 



replacement mg + ag^, where a is an additional quenched chiral parameter. This 

gives 

It is sometimes useful to think of the quenched the limit of the partially quenched 

case as the sea quark masses go to infinity (at fixed valence masses and fixed ml, a, 6'y and 
5^). For disconnected propagators, the resulting decoupling of the sea quarks has the simple 
effect of canceling the "unmixed" terms in the numerator with the terms involving the full 
masses in the denominator. Thus eq. (fTfjl becomes eq. (^Hjl . and eq. (jHIH) becomes eq. (jSH). 
(The aq^ term in eq. pip could have been put in for free in eq. ()30p since it is irrelevant in 
the ml —>■ oo limit.) 



III. ONE LOOP PION MASS FOR 4+4+4 DYNAMICAL FLAVORS 

We can now calculate the 1-loop Goldstone pion self energy. We shall use the term "pion" 
to refer to a generic flavor-nonsinglet meson here, so it can refer to the kaons as well (and 
also what we will shortly call a meson) . As in Ref . Q , all the contributing diagrams are 
tadpoles, as shown in Fig. Q, coming from each of the terms in eq. ()17p . We can break up 
the self energy (defined to be minus the sum of self energy diagrams) as 

^(P') = + ^'''%P')] , (32) 

where "con" and ^^disc" are short for connected and disconnected, respectively. The main 
difference here from Ref. |3] is that the disconnected piece (for m„ ^ m^) now receives 



contributions from all the terms in the Lagrangian, not just t 



le mass term. Also, note that 



we have factored out l/lGvr^/^, and not l/967r^/^ as in Ref. 

The terms "connected" and "disconnected" refer to the internal loop at the quark level. 
In other words, a disconnected diagram will have either an internal disconnected propagator 
(Figs, infg)-(j)) or a disconnected vertex (Fig. Efe)), or both (Fig. Eff)). The disconnected 
propagators correspond to one or more insertions of a two-point 6' or ml vertex {i.e., V in 
eqs. fl77|) . (PH|) and (jHU)) ): while the disconnected vertices are generated by the W term in 
the potential, as we will see below. 

We will explicitly perform the partially quenched calculation. Here the quenched valence 
quarks (call them x and y) will in general have different masses from the sea quarks u, d and 
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s. We will still refer to this as a "4+4+4" partially quenched theory, based on its dynamical 
quark content. The chiral Lagrangian needed has 5 flavors, but 2 flavors are dropped, by fiat, 
from loops. From this 5-flavor, partially quenched theory, we can find equivalent 3-fiavor 
(what we call "full" theory) results by setting the valence quark masses equal to various sea 
quark masses. 

The valence quarks x and y form new mesons in our theory, which we name as follows: 

X = XX Y = yy 

= xy = yx . (33) 

We will not give individual names to the mesons formed from various valence-sea combina- 
tions such as xu, but just refer to them generically by "Q." A check on our flnal calculation 
here is that the 1-loop correction to m'i+ should be proportional to rux + rriy (and hence 
m? . itself), due to the separate Ua{^) symmetries for x and y and the x ^ y interchange 
symmetry. 

Since the mass, kinetic, and U terms are composed entirely of single traces, the relevant 
4- meson vertices that they generate are all of the form of Figs.|3fa) and (b). (This is because 
"touching" flavor indices must be the same in a single trace.) In Fig.|3fb) the vertical meson 
lines must join to make the internal loop; however, they can only join with a disconnected 
propagator because they have different flavors. Thus connected contributions from mass, 
kinetic, and U all involve the vertex of Fig. EJ^a), and produce diagrams of the form of 
Fig-Efa). Disconnected contributions can come from Figs. El^a) and (b). 

The U ' terms, on the other hand, involve two traces, and therefore generate disconnected 
vertices, which in principle can be of the form of either Fig. |3fc) or (d). However it is not 
hard to show, from the explicit taste structure of W, that only vertices with odd number of 
mesons coming from each trace contribute when two of the mesons are Goldstone particles 
(pseudoscalar taste). The W vertices then separate into two disconnected pieces, one with a 
single meson and the other with three. Thus U ' vertices must be of the form of Fig. Efd) , and 
not (c). This in turn implies that the U ' self energy diagrams have the disconnected structure 
of Fig. Efe) or (f) only, where (e) uses a connected propagator and (f), a disconnected one. 

Combining the connected contributions, we flnd: 

"-(f) = E / $ [P' + + ("V + <■) + "'^f^")] ?T^- (34) 

1^ n •' H ' Ob 



Q,B 
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As before, B takes on the taste values {5, /i5, fiu {fi < u), fi, I}, and Q runs over all 
meson flavors with one valence quark {x or y) and one sea quark {u, d, or s). Which 
mesons contribute is clear from Fig. |21^a).^° The first two terms in eq. (j34p come from the 
kinetic energy: one from the derivatives acting on the external legs and the other from the 
derivatives acting on the internal loop. The last two terms are from the mass term and 
U respectively. We have used the fact that /i(ma; + my) = + to rewrite the mass-term 
contribution. 

The one- loop mass renormalization is just the self energy with the external momentum 
evaluated at —m? + . Making this substitution and noting from eq. (fTSl) that = 
+ a^A(^5), the term inside the square brackets becomes + ^q^, which cancels the 
denominator from the propagator. Thus, no chiral logarithms arise from these terms. This 
corresponds to the fact that all diagrams of the form of Fig. Of a) cancel in the standard 
continuum chiral logarithm calculation (See Ref. Q for more discussion.) 

For the disconnected contributions, it will be convenient to divide up a'^^'^'^ further, ac- 
cording to (1) whether the particle in the loop is a vector, axial-vector, or singlet in taste, 
and (2) the type of diagram that generates the term. We thus have 

disc disc , disc , disc /'Ok\ 

with 

(jy'^'^ = av,gh + C''V,ij + '^V,e + CTyj 

^(hsc ^ ^^^^^ _|_ _|_ QT^^^ _|_ ^^^^ 



a 



disc 



(^I,gh + CTl^ij . (36) 



Here, the labels gh, ij, e, and / refer to the diagrams in Fig.Elthat generate the contribution. 
As discussed above, the gh and ij contributions come from kinetic energy, mass, or U vertices, 
with a disconnected propagator. The e and / contributions have a U ' vertex and a connected 
or disconnected propagator, respectively. As is easily seen from the form eq. ()14|1 . U ' vertices 
that have two Goldstone mesons on the external lines must have only taste-vector or axial 
mesons on the loop. Therefore, af'^'^ gets contributions from Fig. |ni^g)-(j) only, as it does in 
the continuum. 



^'^ If we were considering a full n — h flavor theory where x and y were unquenched, then the quark loop in 
Fig. OJa) could also be an x or y, and the sum over Q would include the mesons X, Y and P (eq. 
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We focus first on the taste- vector contributions, cxy^e uses tlie vertex Fig. HJ^d) with i = x 
(or i=y in its x y variant) in order to have a connected propagator. We find, therefore, 



where all but the overall coefficient follows immediately from the form of the diagram. We 
have already included the factor of 4 for the four degenerate taste-vector mesons, and will 
continue to do so below. 

avj again uses the vertex Fig. 111(d), but now i must be one of the sea quarks, since a 
virtual quark loop is involved. The propagator is the disconnected taste-vector propagator, 
V^, eq. (I2ZI)- We have 

2 
3 



^y,/ = -^a%/^ 1^ {-DU + V^) . (38) 

M=U,D,S 

Note that both cry g and ayj have explicit factors of 5'y. These come from the 4- meson vertex, 
generated by W . There are additional implicit factors of 5y in the disconnected propagator 
T)^ in (Tvj- It is not immediately obvious that this same linear combination of C2V and C^v 
(see eq. (|2Uj)) must occur in both the 2-meson and 4- mesons vertices. However, we will see 
below that it is necessary for the cancellations that allow the mass renormalization to 

be proportional to m? as required by axial symmetry. 

^5 

(Tv^gh is generated by vertices of type Fig. ID^a), with i = y {01 i=x in its y — x variant). 
The result is 



71 /f TT r\ a 



+ (p^ + q^ + ml+ + m\ + a^Ay j V 



V 

XX 



V 

YY 



The p^ + q^ terms come from the kinetic energy vertex; m 



m 



Xs 



(39) 

and ^y-g, from the mass 



vertex; and the Ay terms, from U. Putting 
from eq. (jlSp . this simplifies to: 

1 



— m^+, and using nix^ = + a^Ay, 



Wy 



(40) 



Finally, we have (Tv,ij- This contribution uses the vertex Fig. 111(b) and, clearly, an X-Y 
disconnected propagator. 

2 r 



(^V,ijiP 



2 I 2 

p + q 



XY 



2m\+ + Ay 1 



XY 



(41) 
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The sum of all the contributions to ay^'^ can be simplified with an identity derived by 
combining eqs. and (|^: 



M=U,D,S 



V 

XM 



1 



+ m 



XX 



(42) 



Using this and the X Y version we can add eq. (|HHj) to eqs. P7j) and The trivial 

identity (from eq. (EH)) 

{q' + m\^)Vl^ = {q' + ml^)Vly (43) 

and the fact that + my = 2m? + = 2m\,+ + 2a^Ay, can then be used to combine the 
result with eq. (PT|l to give simply 



a 



disc/ 



-m 



2mi 



2 ^XY ■ 



(44) 



Note that the result is proportional to mp+, as expected. The corresponding expression for 
^(ksc obtained hy V A. 

For af***^, we just have contributions from Figs. Etg)-(j). These contributions are very 
similar to the corresponding ones for cxy*'^. We have: 

74, 



+ (p^ + q^ + m\+ + + a^A/ j P^y 



o-i,gh[-mp^ 



12 J n 

A 

6 J TT^ 
1 /"A 



[(g2 + rn^J + (g2 + ^2^J 



(45) 



+ q^ - mi,+ + a^Aj 



XY 



q^ - 2ml+ + a^Ai 



V 



XY 5 



6 J TT^ 

The V I version of eq. (jlHj) allows us to combine eqs. and pUj) . yielding 



(46) 



disc/ 



m 



-m. 



(47) 



Again, the result is proportional to mp+. 



Collecting eqs. pi]) and according to eq. gives 



-^/ jmY+^Viy-V^^y) . 



(48) 
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Since a™" is just a quarticly divergent constant, the above result contains all the 1-loop 
chiral logarithms in the mass renormalization. 

The result in eq. (I48|) is rather implicit. To express the chiral logarithms more concretely, 
we would need three further steps: 

(1) find the explicit expressions for the eigenvalues of the full mass matrices in the 
denominators of the V {e.g., m^o , mf^^ and m^, in V^). 

(2) take the oo limit in the taste-singlet term. 

(3) write the disconnected propagators T> as sums of simple poles and perform the 
integrals over q. 

Steps (1) and (2) are slightly different in the 1 + 1 + 1 case of interest than in the present 
4+4+4 case, so we postpone them until later. On the other hand, step (3) can be done quite 
generally, so we present it here. 

The integrands in eq. ()48p are of the form 

Jl-'^l i{m}M) ^ "r^^^l^l , (49) 

where {m} and {fi} are the sets of masses {mi, m2, . . . , m„} and {/ii, /i2, . . . , /ifc}. respec- 
tively. As long as there are no mass degeneracies in the denominator, and n > k (which is 
true here even after the ml oo limit), I[n, k] can be written as the sum of simple poles 
times their residues: 

j=i y ^ ''l 

where 

4-(0.};M)^ g^g::| . (51) 

Equation ()50|) just follows from the fact that an analytic function is determined by its poles 
and behavior at infinity; it is known as "Lagrange's formula" in complex analysis [iG^ . 
The integrals of the simple poles can now be done using 

^1 = / (tX 2^ 2 ~- TTT^^K) ' (52) 
J (27r)* + lovr^ 



where 



2 

(m^) = In — — - [infinite volume] , (53) 
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with A the chiral scale. We use the arrow in eq. ()52p and later to indicate that we are only 
keeping the chiral logarithm terms. If the system is in a finite (but large) spatial volume 
L^, we only have to modify eq. 

In — + 6i{mL) j [finite spatial volume] , (54) 

where |3| 

with Ki the Bessel function of imaginary argument. 

With the above, we can write a general integral of the form in eq. P8|) as 

/74 TL 
^J["''=l {{myM) - ii^yM) ^K) • (56) 

i=i 

We make one final comment on the 4+4+4 calculation before going on to the 1+1+1 case. 
In Ref . , certain chiral logarithm terms were claimed to come from pure valence diagrams, 
with connected propagators; while in the current calculation, all such terms cancel. What 
is the reason for the discrepancy? As discussed above, the problem in [3| was the incorrect 
treatment of fiavor indices. Because of this, it was not realized that there is a difference 
between a propagator of a fiavor- neutral Goldstone pion, such as f/5, and that of the fiavor- 
nonsinglet vr^ (or between their partially quenched counterparts, and P^)- An explicit 
computation in the current framework shows that the connected, valence terms found in \l\ 
do in fact exist, but only for a fiavor- neutral propagator. Such terms arise identically in the 
Xs-Ys, X5-X5, and Y^-Y^ propagators, but do not appear in the P^-P^ propagator. This 
proves that they come from Fig. Efd). The needed vertex is Fig. EJb) (after relabeling), 
which is generated by kinetic, mass and lA terms. The claim in [7| that connected, valence 
terms come from Fig. Efc) with vertex Fig. |^c) is incorrect. Indeed, it was argued above 
that the fiavor structure of the terms in our Lagrangian forbids vertex Fig. life), at least 
with two external Goldstone mesons. 

IV. MOVING FROM 4+4+4 TO 1 + 1 + 1 DYNAMICAL FLAVORS 

To make the 4+4 + 4 result, eq. ()48j) . into a 1 + 1 + 1 result, we simply must divide by a 
factor of 4 for every sea quark loop. The contributing diagrams are Figs. El (e)-(j). There 
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can be either taste-vector, taste-axial vector or taste-singlet mesons on the internal lines of 
these diagrams, and we can treat all these cases at the same time simply by defining 

{a?5y, taste- vector; 
a^5^, taste-axial; (57) 
4mQ/3, taste-singlet. 

Diagrams (e), (g), and (i) have no sea quark loops and a single factor of 5' (in (e) this comes 
from the 4-meson W vertex). Diagrams (f), (h), and (j) have one additional factor of 5' 
for each sea quark loop. Therefore, dividing by 4 for every sea quark loop is the same as 
dividing every factor of 5', except the first, by 4. For a general function f{6') which vanishes 
linearly as 5' 0, we can make this adjustment simply by the replacement f{5') —>■ 4/ (j). 
Alternatively, we can see from eqs. P8|l . pTjl and (j30|) that the first factor of 6' comes from 
the explicit 6' in front of , T>^ or V^; while higher order terms in 5' are implicit in the 
values of the "full" masses in the denominators relative to the "unmixed" masses in the 
numerators. Therefore, to go from 4+4+4 to 1 + 1 + 1, we leave the explicit 5' factors alone 
but merely let 5' — > 5'/4 before diagonalizing the full mass matrix. 

The full mass matrices to be diagonalized follow from the fiavor-neutral mixing term in 
£, written down following eqs. pO|l and (j2T|) . After 5' — 5'/4, these have the form 

/ml + 5' /A 5' /A 5' /A \ 

574 ml + 5' /A 5' /A (58) 
\ 6' /A 6' /A ml + 6' /A J 

Here the masses mfj, mj^, m| have an implicit taste label (V", A, or J) depending on 
which case we are considering. The explicit expressions for the eigenvalues of eq. ()58|1 are 
complicated and not illuminating in general. The solutions in the 2 + 1 (m^ = m^) case, 
however, have simple forms, and that is the case of greatest current interest. In the taste- 
vector channel, we have, for the 2 + 1 case. 



^I'y = \ (<v + "^Iv + \(^''^'v + Z] ■ (59) 



I 2 2 A2 (^'^^'v I 2 2 \ , %a^^'vf 

\^Sv - Kv) - i'^Sv - Ky) + 



16 



22 



The taste-axial case just requires V ^ A. In the taste-singlet case, S' = 4771^/3, and mg will 
be taken to infinity, so only the large-mo expressions are needed. We have (again for 2 + 1): 

mlo = ml^ = ml^ 

< = ^ + ^ (60) 

mj, = ml , 

where we have neglected corrections that are 0{l/mD compared to the terms kept. 

Finally, we can give the result for the chiral logs in the Goldstone pion self energy. For 
the moment we stay with the partially quenched expression and also assume no degeneracies 
among the valence and sea quark masses. In the 1 + 1 + 1 case we obtain from eq. (jlHj) with 
eqs. (EZ|), (EOl) and (jSl: 

2 

Vl+l+l/' ^2 \ ^ -^5 / 0^2e/ d[5,3] fji2 



jv JA 



2 



+iE4f^K)). (61) 



3 



where we have used eq. (j^ . and R^' and liw?) are given by eqs. and or jy 



runs over {Xy, Yy , iiy , rjv , rj'y} and 



(For Rff\ just let V ^ A). Similarly, ji runs over {X/, Yj, vr^, rji} (the rj'j has decoupled in 
the m-Q — oo limit), and 

= i?^' ({m^,, m^,, m.o, {m^,, m^,, . (63) 

In each taste channel, the values of m^o, and m^, in eq. (jUTjl are just the eigenvalues of 
the corresponding version of eq. 

The 2 + 1 (m^ = m^) case is very similar, but because m^o = "tTT-u = ""^d, there is a 
cancellation in eqs. (jTfj) and (jHUj) between + m^o in the denominators and, say, + 
in the numerators. Assuming no other degeneracies, we have 

2 



ml 
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Here jv runs over {Xy, Yv,r]v,riy} and 

= ^ {{m^.,my^,m,,,m,,J- W,,m.,}) . (65) 
Again, let V ^ A for Rff^- The index ji runs over {X/, Yj,rii}, and 

^ff^ = ^j/^^ m^,}; {muj, msj}) . (66) 

In this case, the values of "^^o, and m^, are given by eqs. and (jUUj) . 

Cases of interest with further degeneracies (such as a "full" 2 + 1 pion with m^. = rriy = 
rriu = ma) can be obtained by carefully taking limits in eq. (j64j) . We will write down some 
of these cases explicitly in Sec. IVll where we also include the analytic contributions. 



V. QUENCHED CASE 

Since we can think of the quenched theory as the limit of the partially quenched theory 
as the sea quark masses go to infinity, all the manipulations that led to eq. (j^Hj) will go 
through unscathed in the quenched case. We can therefore simply replace the disconnected 
propagators in eq. (PHj) with their quenched versions, eqs. and (jHH). Using the same 
notation as in eqs. (jUT|) and (jMj) . we have (assuming ^ m^): 

2 

iv 3A 

h 

Here jv runs over {Xy-, Yv}'-, similarly for ja and j/. For the a-dependent terms, we have 
used the integral 

^= - 1 wr-fh - - / - - • 

where Xi is defined in eq. fjH^ . 

Because the quenched residues here are particularly simple, it is useful to write out the 
result more explicitly: 



2 



j-i quench ^ 



m 



2 ^ , "^^^ 



2e/ ^-y^xA) ~ "-y^YA, 



- °^\v ^ ^Ya - ^Xa 



2 {ml — am\^)l{m\^) — [ml — amY^)£{mY^)'^ 



3 — 



(69) 
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VI. FINAL ONE-LOOP RESULTS 



The mass at one loop is given by 

(m;-;°°P)2 = mj,^ + S(-m5,0 • (70) 

The chiral logarithm contributions to T,{—m'i + ) are presented in eqs. (jHH), (jMj) and 
but for complete one- loop expressions we also need the analytic terms. The latter 

are unchanged from Ref. Ql- However, for the analytic coefficients we now prefer to use the 
more standard Lj, rather than the parameters and K'^ employed in 3]- 
In the absence of any degeneracies, we have, in the 1 + 1 + 1 case, 

l-loop, 1+1+1\2 



+ TUy) 



2 

iv Ja 



+ ^ (2^6 - ^4) (m„ + + m,) + a^cj . (71) 

Definitions here are the same as in eq. (jfill) : the chiral logarithm function l{m?) is given 
by eq. ()53|) . or in finite volume, by eq. Recall that ours is a joint expansion in the 

quark masses (generically m) and a?. The analytic terms in (m^~'°°P)^ here are 0{m?) or 
0{ma^)\ 0{a^) terms cannot enter here because the pion mass must vanish in the chiral 
limit. Lattice effects violating continuum rotational invariance cannot show up at this order 
for the Goldstone pion — see Appendix. 

Similarly, for m„ = = rrii (the 2 + 1 case), but with no other degeneracies, we have 



[m^ + mj^) 



JV ]A 

+ 1 E ^(O) + ^ (2^8 - ^5) (m. + m,) 
32/i 

7^ 



+ ^(2L6-L4)(2m; + m,) + a2c|>. (72) 



Definitions here are the same as in eq. (jMj) . 
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The quenched result is 

1— loop, quench\ 2 



{m^ + my) 



167r2/^ 



9 , ^(^xJ 



{ml) 



m 



m 



Ya 



m 



Xa 



2 {ml — am\^)t{m\^) — {ml — amyji^my^) 
16/i 



'Xi 

^Yj - ^Xj 



(2L^-L'5)K. + m,)+a2C"|, 



(73) 



where the primes on Lg, L^, and C indicate that they may have different values than in the 
unquenched cases. Of course, there is no analytic term involving the sea quarks (2Lg — L'^ 
in the quenched case. 

It is useful to write down more explicit versions of the above results in various limits 



pertinent to many simulations. First, with ^ ma, we set m^ 
obtain the "full QCD" charged pion mass in the 1 + 1 + 1 case: 



[m 



l-loop, 1+1+1\2 



( m„ + md) 



167r2/2 



m 



Sv 



m^o 



2 

V'v 



m,, and m,, 



{^o) + 



md to 



- OKo - ml)'^ + (ml - ml){ml - m^, 



+ (. 
16/i 



A] + 



m 



Si 



m\ 



m 



VI 



-^(<o) + 

.0 

32/i 



2 2 

mq — m„ „ 



m 



m 



VI' 



VI 



P 



{2Ls - Ls) (m„ + mrf) + — f {21^ - L^) (m„ + + m,) + a'C}. (74) 



/ 



This result is most easily obtained by taking the degenerate mass limits in eq. 



before 



the integral is performed, rather than in eq. (jHJ). The quantities m'^o, fnf^ and m^, are 
eigenvalues of the mass matrix, eq. ()58|) . From eq. ()74j) we can get the charged kaon mass 
simply by interchanging the explicit labels d ^ s and D ^ S. (The neutral labels vr*^, t], 
and 7]' are unaffected.) 

The results for the full pion and kaon in the case of degenerate up and down quark masses 
(both set to mi) are also of interest, as they are needed to fit many simulations. Since the 
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pion and kaon results look quite different, we show them both: 

^l-loop, 2+1^2 



2mi 



1 



1 



167r2/2 



1 1 
m% — mt 



rjy Vv nVr Vv 



m\ - ml, 



2 ■ 
V'v' 



+ 



V 



+ ^ (2L8 - L,) {2mi) + ^ (2L6 - L,) (2m, + m,) + a'C 



(75) 



^5 

{mi + ms) 



^ 1 



^«)-^Kj) 



V ^ A 



+^^Kj) + ^ (2L8 - Ls) (m, + m,) 



32/i 

'7^ 



(2L6 - L4) (2m/ + m,) + a^C L 



(76) 



Again, the relevant limits are most easily taken before the integrals are performed. The 7r°, 
7] and rj' masses in this case are given explicitly in eqs. and (jUn|) : we have made heavy 
use of these explicit forms to simplify the chiral logarithm terms in the vr mass. 

The last case we will look at is the quenched pion mass correction in the limit of degenerate 
valence masses (m^, = m^). Here we get a double pole in the pion self energy. We can either 
carefully take the limit my m^ in eq. (fTSj) . or return to eq. ()48|) with quenched T) terms 
and do the double pole integrals directly. We follow the latter approach. We need the 
following integrals: 



X3 = 

where Ti is given in eq. (f52 

iim"] 



d 



(27r)4 (g2 + m^)2 
(27r)4 (g2 + ^2)2 



dm^ 



Xi — m^Xa 



167r2 
1 



m 



167r2 



(m^) — m^i^m^) 



(m 



in eq. (j53|l : and 



1 1 [infinite volume] , 
with A the chiral scale. In finite spatial volume L^, 



m 



where 



m 



In -— + 1 + 6s{mL) [finite spatial volume] , 



63imL) = 2j2Koi\^rnL) 



(77) 
(78) 

(79) 

(80) 
(81) 
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with Kq the Bessel function of imaginary argument. Note that the +1 term in £(m^) was 
omitted in Ref . [3| • That is formally acceptable when we are only keeping chiral logarithms, 
but inconvenient, since then the result from performing the double pole integral is not equal 
to degenerate limit of the chiral logs from the single poles. 
Using the above integrals, we get: 

l"^p+ J r 1 / ~ ~ 9rr?2 _ 



2mx { IGn'^f 



+ - <A<j) (24 - 4) (2m.) + a'C'j, (82) 

Taking the niy limit in eq. (|73p of course gives the same result. To see that the 

finite-size corrections are the same both ways, one needs the identity (7|] 

63{mL) = -6i{mL) - ^6[{mL) . (83) 

Double poles also appear in some other interesting limits of eqs. (17111 and ()72|) . For 
example, the "partially quenched degenerate pion" in either the 2 + 1 case (m. = rriy ^ mi) 
or the 1 + 1 + 1 case (m. = rriy ^ rriu and = rriy ^ met) has double poles. These can be 
dealt with as in the quenched case: either take the limit my in eq. (f7T|) or (f?^ . or 

return to eq. ()48|) and perform the double pole integrals directly. 



VII. REMARKS AND CONCLUSIONS 

The most general result we have is for the n = 3 partially quenched case (1 + 1 + 1) with 
all valence and sea quark masses different, eq. ()71|). Other interesting cases can be obtained 
from eq. (f7T|) by taking appropriate mass limits. The results most relevant to current MILC 
simulations are those with = m^ = mi (the 2+1 case); these and other important limits 
are presented explicitly in Sec. IVII The result in the quenched case is given separately in 
eq. (USD- 

At this point, one can calculate any other desired quantity within this framework. The 
calculation for the pion and kaon decay constants is straightforward; a description is now 



being prepared for publication 



17l |. As in the case here of the one loop pion mass, it is again 



If one chooses to perform the double pole integrals directly, eq. H50() is no longer valid, and a generalization 
of this formula is needed. 
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simpler to examine the partially quenched case, and from there all the necessary results 
can be obtained. The next step will be the incorporation of heavy quarks, so that we can 
examine the effects of staggered discretization errors on heavy-light meson quantities. This 
requires an extension of these ideas to incorporate the heavy quark symmetries within SXPT, 
and is in progress. 

The generalization of the Lee-Sharpe Lagrangian to multiple flavors has shown that two 
additional parameters, 5y and 5^, appear in the one- loop chiral logarithms for the charged 
meson masses. These parameters are not determined at tree level by existing lattice data 
for pion mass splittings, since they contribute only to unmeasured disconnected tree graphs 
for flavor-neutral, taste- nonsinglet, pions. The new parameters are therefore unconstrained 
in current chiral- logarithm fits to lattice results. In contrast, the masses of the charged 
pions of various tastes that appear in our final results are not free parameters in the one- 
loop fits, since they are determined at tree-level by lattice measurements. Using tree-level 
information, a fit of lattice data to eq. (f72|l would have 7 free parameters: /z, /, 2L^ — L5, 
2Lq — L4^, 6y, 6'^ and C}"^ We remark that existence of the parameters 6'y and 6'^ leads to 
the possibility of phase transition before the chiral limit of the staggered theory is reached. 
This possibility is discussed further in the Appendix; it does not appear to be realized in 
practice for the strange quark mass at its physical value. 

Despite the presence of additional parameters, well controlled simultaneous fits to par- 
tially quenched lattice results for /^r, fx, ^l/i^mi) and m\/ {mi+ms) at fixed lattice spacing 
appear possible These should allow for highly accurate extrapolations to physical quark 
mass and then to the continuum, as well as determinations of the Gasser-Leutwyler param- 
eters Li. It can help here to constrain, at least weakly, the new chiral parameters. One easy 
way to do this is to use a vacuum saturation estimate of the matrix elements of the 4-quark 
taste- violating operators calculated in perturbation theory More accurate lattice eval- 
uations of the matrix elements, or perhaps even direct lattice determinations of the 6y and 
5^ by evaluation of disconnected pion propagators, may also be envisioned. 

An alternative approach to the fitting of lattice data is also possible when highly accurate 
data exists at more than one lattice spacing. Here one can extrapolate to the continuum at 



One may choose to absorb C into which wiU have 0{a?) corrections in any case from higher operators in 
the effective continuum action that have the same symmetries as the lowest order terms — see Appendix. 
However, this will change the higher order corrections to eq. (|72|l . 
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fixed quark mass and then fit the resulting "continuum" results to standard XPT forms, i.e., 
without taste violations. This is the approach taken in j^^, and it works well. Because of the 
nonanalytic dependence on the lattice spacing induced by the chiral logarithms coming from 
pions of various tastes, though, there is a residual discretization error left in the data even 
after extrapolation to the continuum. This error would go away if one worked very close to 
the continuum limit, where "very close" here means ka'^Aq^j^ <^ m^/A^, with k a constant 
that depends on the particular staggered action, and A is the chiral scale. For pions light 
enough for XPT to be applicable, however, this condition is very difficult to satisfy without 
further improvement in the staggered action than is currently available. The SXPT formulas 
above will therefore remain crucial, at least in the near term, for determining the systematic 
errors in the results. 
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APPENDIX 

Here we write down the symmetries of the effective continuum action ( "Symanzik action" ) 
of the staggered lattice theory through O(a^), and those of the corresponding chiral theory, 
eq. (fTTj). We also briefly discuss the interesting possibility of a transition of the staggered 
theory to an unusual phase. We follow the notation and reasoning of Ref. jol closely; the 
discussion in this Appendix is not self-contained. 

The symmetries of various terms in the O(a^) Symanzik action are shown in Table HI 
which is a generalization of Table 1 in Ref. Q to the current ra-flavor case. 

The "residual chiral group," U{n)i x U{n)r, which is a symmetry of S^^^^^ and S'^^*-^'', 
is the extension to multiple flavors of the residual f/(l)vcc x U{1)a symmetry of a single 
staggered field. Let ta be the U{n) generators, and let q be the complete (fiavor ® taste ® 
spin) quark field, as in eq. (jH)) but with fiavor indices suppressed. Then the residual chiral 
group is given by: 
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Term in action [Flavor] x Rotation symmetry 



54 (m = 0) [C/(l)vEC X SU{4n)L x SU{An)ji\ x 50(4) 

Si (m ^ 0) [( [7(l)vec X SU{AU, )"] x SO{4) 



oglue 


[?7(l)vEC X SUi4n)L 


X 5t/(4n)R] 


X 5VF4 


^bilin ^ 0) 


[C/(1)VEC X SUi4n)L 


X SU{4n)R] 


X SWi 


^bilin _^ 0) 


[( C/(l)vec X 5C/(4),ec 


)"] X 5W4 




c.FF(A) 
•^6 


[;7(n)f X U{n)r x (r4 


XS0{4))] 


X S0(4) 


oFF{B) 
•^6 


C/(n)^ X U{n)r X (r4 







56 (m = 0) U{n)i x ?7(n), x (r4 >^ 5VF4,diag) 

Seim ^ 0) ( ;7(l)vec )" X r4 Xl 5W^4,diag 



TABLE I: The flavor and rotation symmetries respected by various terms in the effective action. 
Here "flavor" is used generically to include fermion number, true vector flavor, chiral, and taste 
symmetries. Almost all the notation is from Ref. ^6']. The "residual chiral group," U{n)£ x U{n)r 
is defined in the text. We have also added the subscript "vec" for vector (L + R) symmetries, 
and have included overall fermion number, C/(l)vEC5 ^ well individual flavor number symmetries, 
U{l)vec- There is no clear separation of flavor and rotation symmetries in the last three lines. 
For simplicity in the m ^ cases, we assume that all quark masses are nonzero and different for 
different flavors. 



where 91 and 9^ (a = 1, 2, 



exp 


i91K (- 


q exp 


-l9pa 


exp 




q exp 


-i9':.ta 



1 - 75 ® 6 



I + 15® it 



1 + 75 ® 6 



1 - 75 ® 6 



(84) 



v?) are the group parameters. We use the notation I and r, 
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rather than the usual L and R for chiral rotations, because these symmetries combine chiral 
spin with taste. To study the effect of this symmetry on various terms in Table lU consider 
a flavor- singlet, "odd" bilinear, i.e., a bilinear of the form 9(75 ® ^t)"?, where there is an 
implicit sum over flavor, and where {75 ® ^t-, 75 ® ^5} = 0. Then it is clear that any bilinear 
of this type is invariant under the residual chiral symmetry, eq. (|84p . That the S'^^*-^-' and 
5'FF(B) terms in the action are invariant under this symmetry now follows from the fact they 
can be written as sums of products of such bilinears (see discussion preceding eq. (O)- 

Note that, even though the identity matrix in flavor is included among the generators ta 
in eq. (fH^ . the corresponding axial symmetries {01 = —6*") are traceless in flavor ® taste 
because of the presence of ^5. 

For future purposes it is convenient to rewrite eq. (j84j) to show explicitly the action of 
the i and r symmetries on the chiral fields. Define 

1-75 



Ql = 
Qr = 
UiL = exp 
UrL = exp 



2 

1+75 



Ql 
Qr 



-/ I +75 

1 - 75 



(85) 



2 

1-^5 



UgR = exp 
UrR = exp 



^0% 



^5 



2 

1+^5 



(86) 



Then 



Ql 
Qr 
Ql 
Qr 



UiL Ql 
UeR Qr 
UrL Ql 
UrR Qr 



Ql 
Qr 
Ql 
Qr 



QL Ul 
Qr U}r 
QL Ul^ 



QrUI 



rR 



(87) 



One now assumes that the SU{An)L x SU{4n)R approximate symmetry [i.e., the sym- 
metry of 5*4, the 4-dimensional terms in the action, at m = 0) breaks dynamically in the 
usual way down to SU{An)vcc- The kinetic energy term in the effective chiral Lagrangian 
then has the complete SU{An)L x SU{An)R symmetry (realized nonlinearly) . Other terms 
in the Symanzik action are represented by additional terms in the chiral Lagrangian with 
the corresponding symmetries. 

A key insight of Lee and Sharpe is that the chiral representatives of all terms in the 
action that violate the SO {A) rotation symmetry must contain derivatives. For example, 
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the rotationally noninvariant term ® I)D^^q in 5*3''™ has a lowest chiral repre- 

sentive Tr((9^S(9^S''"). Chiral terms that are already 0{a'^) and also have derivatives 
will be higher order than our 0{m,a?) Lagrangian, eq. (jl7|) . Thus only S'4(m = 0), 
5*4(772 7^ 0), and 5*5 contribute to eq. (jl7|) . giving the kinetic energy, mass term, and 
potential V, respectively. The symmetry group of the chiral Lagrangian is therefore sim- 
ply the intersection of the symmetry groups of the relevant three lines from Table IH^^ 
[( f/(l)vcc )" X XI SO{A))] X 5*0(4), although by treating the mass and taste violating 
matrices as spurious, one can work with the full [f/(l)vEc x SU^'in)^ x SU{An)j{\ x 5*0(4) 
group. 

Under the the residual chiral symmetry U{n)^ x U{n)r, eq. (j87|) . the chiral field E trans- 
forms as 

i: T.^U,L^ U}^ , St ^ Urn 

r: E^UrL^Ul^, St ^ f/,^, St f/J^ (88) 

It is straightforward to check that the kinetic energy and potential terms in the Lagrangian, 
eq. (fT7j) . are invariant under this symmetry, which is of course violated by the mass term. 

Note that terms violating continuum rotational invariance can appear for the first time 
at 0{ma ), from the chiral representatives of 5g . However, because the taste of the 
Goldstone pion transforms trivially under lattice rotations (5*W4,diag), rotational violations 
cannot affect it unless four derivatives are present. Thus, for example, the Goldstone pion's 
continuum dispersion relation is violated at (9(m^a^) by a term o^^fj_d'jji^ d'jji^ coming 
from the chiral reprentatives of the noninvariant terms in 5*g''"^ and Sf^'^. 

The 5*0(4) part of the taste symmetry of the lowest order chiral action guarantees that 
the approximate spectral degeneracies found in Ref. joj persist in the 72-flavor theory. This 
symmetry is "accidental" in the sense that it is not obeyed by the full lattice action and will 
be violated at next order. Note that the taste 5*0(4), and in fact the accompanying discrete 
Clifford group appear only once, as an overall taste symmetry affecting all flavors, and 
not as individual groups for each flavor separately. This can be seen from the structure of the 
four-quark operators, eq. It is related to the fact that the symmetry of the underlying 
lattice theory has just a single r4 x 5PV4diag factor, which is generated by single site 



13 



We are ignoring the discrete symmetries of parity and charge conjugation here. 
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translations and lattice rotations. It clearly cannot act on different ffavors separately since 
the gauge fields must also be translated and rotated. We remark further that, if there were 
r4 XI SO (4) taste symmetries for each flavor separately, they would forbid taste-nonsinglet 
hairpins graphs like Fig. ^ 

Lee and Sharpe have discussed the possibility of an unusual "Aoki-phase" for staggered 
fermions that could occur if the mass squared of one of the non-Goldstone pions vanished 
before the chiral limit. However, since the splittings, A{^b) in eq. ()19|) are all positive for 
existing staggered actions, this scenario seems unlikely to be realized in practice. 

The current work suggests another possibility for an unusual phase: from eq. (j^HI), if S'a 
or 6y is negative and sufficiently large in magnitude compared to the Goldstone masses and 
to the corresponding splittings, or Ay (eq. (fTUj) ). then m^^ or m^^ could vanish before 
the chiral limit. This possibility seems to us not as remote as the previous one, because 
chiral logarithm fits [l8| to existing MILC data tend to give a negative value for S'^ that 
is comparable in magnitude to A a- Taking rriu = ma, eq. (j59|l implies that m^^ vanishes 
before the chiral limit (m„ = rrid = 0) is reached if 

< . -4A. ^^3^,^ (89) 

For the s-quark mass at its physical value, m|g ~ (TOOMeV)^. On the "coarse" (a ^ 0.13 
fm) MILC lattices, a^A^ ^ (275MeV)^. This means that the transition could occur with a 
physical strange quark mass only for 5^ < — 1.9 A a, which does not appear to be satisfied by 
the chiral fits. Further, as a decreases, the fit values of 5^ seem to move further from 5^ 
However, the transition appears considerably more likely to be realized in the unphysical 
case where all three quark masses get small. There 5^ ^^.^^ = — 4Ayi/3, which is comparable 
to fit values of 5^. More study of this interesting possibility is clearly warranted. 
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(a) 



u d 

u d 

(b) 

FIG. 1: The two-point mixing vertex (among taste vectors) coming from the new W term, (a) 
corresponds to the chiral theory, (b) shows the corresponding quark level diagram. We also have 
U-S and D-S mixings and diagonal terms {U-U etc.). There are similar vertices among the axial 
tastes (with a^6'y — a^(5^), as well as the singlet tastes (with a^Sy — Auiq/S). 




(a) 




(b) 

FIG. 2: The only diagrams contributing to the flavor-nonsinglet meson self energy, (a) includes 

all contributions where the propagator in the loop contains no two-point vertex insertions, (b) 
subsumes the graphs which have disconnected insertions on the propagator. The cross represents 
one or more insertions of either the m-Q, S'y, or (5^ vertices. 
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FIG. 3: The quark level diagrams that could contribute to the 1-loop meson self energy. Diagrams 
(b) and (c), which require vertices of the form of Fig. El^c), do not occur in the case of interest. 
Similarly, (d) only contributes to flavor-neutral propagators. Note that (f), (h) and (j) correspond 
to (e), (g), and (i), respectively, with iteration of either mQ, d'y, or S'j^ vertices. These diagrams are 
to be taken as including any number of iterations, thus multiple internal quark loops. 




(d) 



FIG. 4: The possible quark level diagrams for 2^2 meson scattering, where one incoming and 
one outgoing particle (shown horizontally) are fixed to be valence mesons, P"*" = xy. The indices 
i and j represent arbitrary quark flavors. There are two additional diagrams (not shown), which 
are like (a) and (d) but have the roles of x and y interchanged. 
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